Evolution of entanglement under echo dynamics 



Tomaz Prosen^, Thomas H. Seligman^ and Marko Znidaric^ 
Physics Department, FMF, University of Ljubljana, Ljubljana, Slovenia 
^ Centro de Ciencias Fisicas, University of Mexico (UN AM), Cuernavaca, Mexico 

(Dated: February 1, 2008) 

Echo dynamics and fidehty are often used to discuss stabihty in quantum information processing 
and quantum chaos. Yet fidehty yields no information about entanglement, the characteristic prop- 
erty of quantum mechanics. We study the evolution of entanglement in echo dynamics. We find 
qualitatively different behavior between integrable and chaotic systems on one hand and between 
random and coherent initial states for integrable systems on the other. For the latter the evolution 
of entanglement is given by a classical time scale. Analytic results are illustrated numerically in a 
Jaynes Cummings model. 
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More than a century ago J. Loschmidt, in his discus- 
sions with L. Bolzmann illustrating irreversibility in a 
gas, suggested to invert the velocity of each atom indi- 
vidually in order to revert to the initial situation. Re- 
cently Loschmidt echoes have been of great interest in 
the control of quantum information processing 1]. As 
entanglement is the essential property of quantum me- 
chanics, in the present paper we analyze the Loschmidt 
echo appearing in the entanglement of two quantum sys- 
tems. As a measure of entanglement we use purity 2] and 
show analytically as well as numerically that for classi- 
cally integrable systems the purity decays as 1 — cjnt^, 
whereas for a classically chaotic system the decay after 
the Zeno time scale is described by 1 — Cch^- For a coher- 
ent state the constant cin is independent of h. For chaotic 
systems (or random initial states in integrable systems) 
this constant on the other hand is proportional to 
thus defining quite different time scales. 

Zurek proposed to use the rate of decoherence as 
a characteristic of chaos in quantum mechanics and this 
is occasionally re-interpreted in terms of entanglement 
Is] between different parts of a closed quantum system. 
Such studies were limited to forward time evolution. Yet 
the insensitivity of quantum mechanics to small changes 
in initial conditions has been a basic difficulty in the 
introduction of the concepts of chaos to this field, and 
the idea to use sensitivity to perturbations in the time 
evolution 4] has emerged recently as one of the tools 
of choice to overcome this difficulty 5, 6]. Specifically 
all these authors use fidelity, i.e. the correlation func- 
tion between a quantum state evolving under the action 
of two Hamiltonians differing by a perturbation, which 
is equivalent to the autocorrelation function under echo 
dynamics. Since fidelity measures irreversibility of a full 
quantum state under the echo, it is also desirable to un- 
dertand irreversibility of a less restrictive quantity like 
entanglement. A recent study in spin chains 7], more 
related to quantum computing, revealed no qualitative 
difference between fidelity and the evolution of entan- 
glement under echo dynamics as measured conveniently 
by purity there denoted as purity -fidelity. This system 
allows no classical analogue, and by consequence no co- 



herent states. Yet these we shall show to be essential to 
recover results analogue to the ones of Zurek and Nemes 
0,0. Based on the idea that a partial trace simulates de- 
coherence our results lead to the conjecture that Zureks 
result holds exclusively for coherent states in the central 
system, i.e. that we may not expect faster decoherence 
for chaotic central systems than for integrable ones, if 
we use random states that are more relevant to quantum 
information. 

To implement this idea we have to consider systems 
with at least two degrees of freedom to allow for entan- 
glement and a well defined classical limit ^ 0. The 
(unperturbed) Hamiltonian will contain a parameter per- 
mitting a transition from order to chaos, and will typi- 
cally couple the two degrees of freedom. A perturbation 
is then defined to obtain a second similar Hamiltonian. 
We give general results for entanglement under echo- 
dynamics starting from an initial product (dis-entangled) 
state. To illustrate our results we use the Jaynes Cum- 
mings (JC) model 8]. The usual co-rotating (integrable) 
version of this model has great practical importance in 
atomic physics and illustrates the h independent evolu- 
tion of entanglement for coherent states. For the argu- 
ments involving the chaotic dynamics we include counter- 
rotating terms 3, 9] to construct a toy model that allows 
for chaos. Even this model may not be entirely unrealis- 
tic for atoms in a Paul trap in a driven field, as standard 
papers seek conditions where this term is small lOj. 

For general considerations and analytic calculations 
techniques of linear response developed originally for the 
evaluation of fidelity 6] are extended to calculate purity- 
fidelity in terms of time correlation functions of the per- 
turbation. In the case of coherent states we could carry 
the evaluation of linear response one step further using it 
in a semi-classical framework that relates the decay rates 
directly to the stability matrix of the orbit along which 
the packet evolves. 

We consider the unitary time evolution given by the 
echo operator Ms(t) = Ul(t)U(t). Here U(t) is generated 
by some unperturbed Hamiltonian H diS U {t) = e~*^^/^ 
and similarly Us{t) = e~*(^+^^)^/^, where V is the per- 
turbation with strength S. It is useful to rewrite the echo 
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operator as time-ordered product ^ in the interaction 
picture 



Ms{t) =f exp (iS(t)(5/n), 



(1) 



where := V(r)dr with V{t) := /7t(t)l/[/(t). This 
operator shah act on a composite system with the Hilbert 
space H = Til ^ '^2, consisting of two factors with di- 
mensions Ni and A/'2, which we may look upon as a 
"central system" and an "environment". We are inter- 
ested only in information about the subsystem 1 which 
is contained in a reduced density matrix pi{t) := tr2p(t), 
p{t) = We shall study the purity-fidelity 7] 

Fp(t) := tri[pi(t)]^ as a measure of factorizability of a 
joint state = Ms(t)\ilj{{))) . We choose this quantity 

rather than some entropy because of its simple analytic 
dependence on pi{t). Here the partial traces with indices 
1 and 2 are taken in the corresponding factor spaces and 
the reduced density matrix is acting on the first factor 
space. We always assume that we start with a factor ized 
state at t = 0, i.e. Fp{0) = 1. For comparison we shall 
also use the fidelity = |(V^(0)|M5(t)|V^(0))p. 

Expanding the echo operator in ^, we get jSf 



\F{t)f = i-s^h-^c{t)^- 

C{t) := - imf^ 



(2) 



Here (.) denotes an expectation in the product initial 
state I '0(0)) = |1, 1) with the abbreviation |z, i^) := \i)i 
\u)2' The same techniques yield for purity-fidelity 

Fp{t) = l-2S^h-^{C{t)-D{t)}^--- , (3) 
D{t) := ^Kl,HS(t)|l,l)p + ^Kz,l|S(t)|l,l)p 



For both series to converge it is sufficient to use a bounded 
perturbation operator V, but we expect the linear re- 
sponse formula to be a good approximation for a much 
wider class of perturbations. The somewhat unusual cor- 
relation function D{t) contains only off-diagonal matrix 
elements of the operator and determines the differ- 
ence between Fp{t) and \F{t)\^. From expansions (|2I3|) 
we can see that the decay is determined by time corre- 
lation functions of the perturbation. The stronger the 
decay of correlation functions {ilj\V{t)V{t')\ilj) as \t — t'\ 
grows, the slower is the increase of C{t) and the slower 
is the decay of F{t) and Fp{t). 

We limit our discussion to systems which have a classi- 
cal limit. For such systems chaos typically implies decay 
of the time correlation functions of the perturbation ob- 
servable {i.e. mixing), while regular motion implies non- 
ergodic behavior. Fidelity decay for both situations is 
discussed in Ref. 0. Under rather general assumptions 
one finds exponential decay for chaotic dynamics 



\F{t)\^ = exp i-t/Tem). Tem = (2(j)-^ HH' 



(4) 



where a diffusion coefficient a := lim^^oo C{t)/{2t) is in- 
dependent of the initial state |'0(O)) [for sufficiently long 



times, typically t ^ logl//z]. In classically re^i^/ar situa- 
tion the fidelity exhibits a quadratic decay in the leading 
order in 5 even for long times, since C{t) ct'^^ where 
c depends on the structure of the initial state. For a co- 
herent initial state we find a Gaussian decay of fidelity 



|F(t)|2=exp(-(t/r„e)'), 



(5) 



with c (X h 6\. It is worth to stress that in the regime of 
linear response (small S) formulae (|4I5|) agree with Eqs. 
(|2I3|) from which the time scales rem,Tne are obtained. 

As for purity-fidelity of chaotic systems, one may argue 
that T^{t) should look like a random matrix so the term 
D{t) should be small compared to C{t) in 0, namely 
D(t)/C{t) ^ 1/Ni^l/N2 because of the smaller number 
of terms involved in the sums. Thus, if both dimensions 
Ni^2 grow as ^ one expects in the asymptotic regime 
that Fp(t) = \F{t)\'^ = exp(— 2t/rem), and does not sig- 
nificantly depend on the initial state. In particular this 
also holds for coherent initial states. Using similar ar- 
guments for regular dynamics but a random initial state 
one agam sees that Fp{t) follows 1^(^)1"^ closely 7]. 

Yet for coherent states and regular classical dynamics 
this is not the case because the term D{t) is not negligi- 
ble. We show that the difference C{t) — D(t) cancels in 
the leading order in fi, i.e. C(t)—D(t) ~ ^ meaning that 
Fp(t) as compared to \F{t)\^ decays on a qualitatively 
longer, fi-independent time scale r^^ = K/S 

^ /T/ ^ne • 

This will be the main result of the present paper. In 
order to establish this we consider the evolution of a 
Gaussian wave-packet along a stable orbit Zt = {xt^pt) 
as {x\ilj{t)) Cexp{^[{x-Xt)-At{x-Xt)^Pfx]), where 
the block-form of the complex d x d matrix 



All Ai 
A21 A 



22 



(6) 



corresponds to obvious division ofd = di-\-d2 dimensional 
configuration space into di and ^2 dimensional parts. At 
is a ratio of two pieces of a classical monodromy ma- 
trix Til SO it is fi-independent. The purity of a reduced 
wave-packet p(xi, x'l) = / dx2{xi, X2\il^{t)){ilj{t)\x[, X2) is 
Fp = J dxidx[\p{xi,x[)\'^ = 1 if A12 = A21 = while in 
general we find fi— independent expression 



Fp = (det ImA) 



Im All 

i J* 
2^21 



1 4* 

2^12 

Im A22 

-M12 






-M21 
Im All 

i A* 
2^21 



--2A12 



Im A22 



-1/2 



For classical echo dynamics, the covariance matrix A = 
At is given by a linear stability analysis as A = + t6B 
for some matrix where Ao,i2 = ^0,21 = 0. Then 
purity-fidelity is fi-independent and can be evaluated in 
the leading orders as Fp(t) = 1 — (tJ/ilT)^ + — 

We thus reach the following interesting conclusion: 
Both fidelity and purity-fidelity decay quadratically in 
integrable situations, while they decay linearly in chaotic 



3 



ones, once we are beyond the Zeno time scale. Yet there 
is a very relevant difference in time scales themselves, 
if we discuss the purity of coherent rather than random 
states. For integrable systems purity- fidelity decays on 
an fi-independent scale. This leads to situations with 
very stable purity-fidelity, while the same perturbation 
generates decay of the fidelity of the coherent state as 
well as the decay of the purity-fidelity of a random state 
on much shorter time scales, dictated by the value of 
h. Note though that for sufficiently small perturbations 
at fixed h the quadratic decline of purity-fidelity always 
prevails. 

To illustrate these results we use the JC Hamiltonian 
including co- and counter-rotating terms for the chaotic 
case as 

H = huja^a + heJ^ + — ^(GaJ+ + G'aJ- + h.c.) (7) 
V 2J 

with standard boson operators a, ^ [a, a^] = l^ and stan- 
dard SU(2) generators J±^Jz. We choose h = 1/ J ensur- 
ing that the classical limit is reached for J ^ oo while the 
angular momentum hJ = 1 is fixed. If either G = or 
= the model is integrable with an additional invari- 
ant being the difference or the sum of quanta for the spin 
and the oscillator. In all calculations we used coherent 
initial states for the product system, i.e. direct product of 
coherent states of the oscillator, \a)2 = e"^ ~" ^10)2, and 
of the spin [SU(2)], \0,6h = (1 + rr*)-'^ exp (r J_)| J, J) 
with r = e^^tan(l9/2) [H. 

For our numerics we fix J = 4 and choose initial posi- 
tion of SU(2) coherent state at (^,0) = (1,1) and for the 
oscillator at a = 1.15. The parameters in JC Hamil- 
tonian are: (a) in chaotic regime co = e = 0.3 and 
G = G' = 1, (b) in integrable regime a; = e = 0.3 and 
G = 1, G' = 0. The corresponding classical Poincare sec- 
tion shows a single practically ergodic component in the 
chaotic case (a) [at energy E = 1.0 determined from the 
initial condition], whereas integrable case (b) [E = 0.63] 
shows a generic family of invariant tori. The perturba- 
tion is realized by varying the parameter e in JC Hamil- 
tonian {TJ, also known as dephasing, so the (bounded) 
perturbation generator is F = hJ^. 

We now show numerical results obtained by diagonal- 
ization in truncated Hilbert spaces. Stability of the cal- 
culation with respect to truncation was checked. Fig. Q 
presents the correlation integrals C{t)/t and D{t)/t for 
chaotic and regular regimes. For chaotic dynamics (case 
a) the correlation integral converges after t ~ 10 to a 
well defined diffusion coefficient a = 0.10 with the D- 
term being of order 1/Ni + I/N2 ^ 1/4. For regular 
dynamics (case b) and t > 10 the correlation integral 
grows as due to a non- vanishing plateau c = 0.046 
in the correlation function. In this case the difference 
G(t) — D{t) is approximately C{t)/J ex fi^, which has 
been checked numerically also for larger J < 24, con- 
firming fi-independent decay of Fp{t). The oscillations 
in these functions are not accounted for by the present 
theory, and are probably particular but interesting prop- 
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FIG. 1: Correlation integrals C{t)/t and D{t)/t for regular 
(top two dotted curves) and chaotic regime (lower two solid 
curves). In both cases the upper curve is for C{t)/t and the 
lower for D{t)/t. The horizontal dashed lines indicate 2a = 
0.20 (upper) and 0.20/4 (lower), whereas the increasing ones 
have the slopes c — 0.046 (upper) and (1 — 0.98/4)c (lower). 
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FIG. 2: Purity-Fidelity Fp(t) (main figure) and squared fi- 
delity |F(t)|^ (inset) in chaotic regime (solid curves) and in 
integrable regime (dotted curves), for S = 0.1. The dashed 
lines indicate the linear and quadratic approximation respec- 
tively. Note the differences in vertical scales. 



erties of the model. Whether they relate to oscillations 
seen in is an open question. 

We first report a calculation with a strong perturba- 
tion ^ = 0.1, that rapidly exceeds the realm of validity of 
linear response, in Fig. [21 where the main figure gives the 
purity-fidelity and the inset the fidelity. For the fidelity 
decay (inset) we find excellent agreement with the expo- 
nential decay in a chaotic regime and a faster Gaus- 
sian decay in a regular regime Q, where the decay rates 
are fixed as above. However, for purity- fidelity we find al- 
ready at t 20, that the decay starts to be influenced by 
the saturation value of Fp{t 00) ^ 1/(2 J + 1). There- 
fore purity-fidelity is higher for the integrable case than 
for the chaotic one not only at short times, as expected, 
but even at large times. This is relevant, because we 
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FIG. 3: (Color) Echo dynamics for weak coupling: S = 0.005. 
Square of the Wigner function for chaotic dynamics (top dia- 
grams) and integrable dynamics as a function of time (bottom 
diagrams) at times corresponding to the axis. Color code: top 
left. Purity- fidelity is shown in the frame on the same time 
scale and for short times in the inset. Red curves give the 
integrable and blue curves the chaotic evolution. Full curves 
show the complete numerics, symbols the evaluation starting 
from the numerical correlation functions of Fig.^and dashed 
curves the linear or quadratic approximation. 

shall next choose a weak perturbation S = 0.005 to avoid 
this problem. We expect and find the crossover after 
a fairly short time. This calculation allows comparison 
with theory as well as an illustration of the evolution of 
the square of the Wigner function , corresponding to the 
reduced density matrix pi (t) for the angular momentum 
states on the sphere using the definition of Near 
the top and bottom of Fig. [31 we see this evolution for 
the chaotic and the integrable Hamiltonian respectively. 
In the center of the figure we plot the purity-fidelity on 
the same time scale as the Wigner functions in the main 
frame and an amplification of short times in the inset. 
We observe detailed agreement of numerics with results 
obtained from the numerical values of the correlation in- 



tegrals (|2I3|) reproducing the oscillatory structure. From 
the same correlation integrals we obtained the coefficients 
for the linear and quadratic decay, which agree well if we 
discard the oscillations. We see a crossing of the two 
curves at t = tp ^ 12 for Fp. This times are larger 
than the Zeno time (^ 1) and indicate the competition 
of the decay rate and the decay shape as expected for a 
non-small value of = 1/4. It is important to remember 
that the integral over the square of the Wigner function 
gives the purity and therefore the fading of the picture 
will be indicative of the purity decay. On the other hand 
the movement of the center is an indication of the rapid 
decay of fidelity (not shown in the figure). 

In this letter we study the linearized behavior of the 
evolution of entanglement under echo dynamics for time- 
scales large compared to those of the quantum Zeno ef- 
fect, but sufficiently short for the expansion to be valid. 
Similarly to the behavior of fidelity the decay of purity- 
fidelity is typically quadratic for non- mixing systems, and 
linear for mixing ones, the first situation arising for inte- 
grable systems and the second for chaotic ones. An inter- 
esting particular, but relevant case appears if we consider 
coherent initial states and integrable classical dynamics. 
In this case we have shown that purity-fidelity, still hav- 
ing a quadratic decay, can be computed classically in the 
leading order which is fi-independent, so the time scale 
for purity-fidelity decay of a coherent state is longer by 
a factor proportional to h~^^'^ than the corresponding 
one for a random state. Coherent states in integrable 
systems are thus particularly long lived for semi-classical 
echo situations. On the other hand, for chaotic classical 
dynamics and coherent initial states we find that purity- 
fidelity is the same as for random states, and its decay 
will be slower than for either random or coherent ini- 
tial states and integrable dynamics provided that time is 
sufficiently long or perturbation 5 is sufficiently small. 
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